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EQUATIONS FOR THE INDUCED VELOCITIES NEAR A LIFTING ROTOR 
WITH NONUNIFORM AZIMUTHWISE VORTICITY DISTRIBUTION 
By Harry H. Heyson 

SUMMARY 


Equations, which can be integrated on high-speed computing machines, 
are developed for all three components of induced velocity at an arbi- 
trary point near the rotor and for an arbitrary harmonic variation of 
vorticity. Sample calculations for vorticity which varies as the sine 
of the azimuth angle indicate that the normal component of induced 
velocity is, in this case, uniform along either side of the lateral 
axis . 


INTRODUCTION 


The mutual interference between the rotors, wing, and tail of a 
helicopter or convertiplane may have a large effect upon the overall 
performance and stability of the aircraft. Since mutual interference 
may entail significant performance penalties, it is necessary to con- 
sider such effects in the initial stages of design. The accurate 
assessment of mutual interference, however, requires a knowledge of the 
flow field of a lifting rotor. The most complete study of the induced 
velocities near a lifting rotor is that of reference 1, which shows that 
the flow may be calculated with reasonable accuracy provided that a 
representative distribution of disk loading is assumed. The theory of 
reference 1, however, is based upon existing induced velocity calcula- 
tions. Since these calculations consider only the case of uniform 
azimuthwise vorticity, the theory, as given in reference 1, can only 
be used to represent disk loadings for which the circulation distribu- 
tion is circularly symmetrical. As a consequence, the flow field, 
according to reference 1, is symmetrical about the longitudinal plane 
of symmetry of the rotor. The measured flow field of reference 1, 
however, shows significant differences in the flow on opposite sides 
of the rotor. Thus, it is desirable to examine the possibility of 
extending the available theory to include the varying components of 
vorticity which are known to exist in the wake of the rotor. 
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The calculation of the flow field is greatly complicated by the 
inclusion of terms representing azimuthwise variations in circulation. 

The only investigation which has considered these effects is that 
described in reference 2, which examined the induced velocity distribu- 
tion on the lateral axis of the rotor for sin t circulation. The 
actual calculation was based upon a crude approximation to the wake that 
is believed to correspond to this distribution of circulation; further- 
more, it is not possible to extend the analysis of reference 2 to points 
other than those which lie on the lateral axis. Thus, it is necessary 
to consider the problem from a different viewpoint if completely general 
results are to be obtained. L 
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The present paper develops equations for all three components of 9 

induced velocity at an arbitrary point near the rotor and for an arbi- 7 

trary harmonic of the azimuthwise distribution of circulation. This 
derivation is accomplished without recourse to the approximations of 
reference 2. Unfortunately, the final integration (with respect to \|r) 
is not possible in closed form; however, numerical results may be 
obtained by the use of modern high-speed computing machines. As an 
example, the distribution of the normal component of induced velocity on 
the lateral axis is computed for sin i|r circulation. These calculated 
velocities are compared with those of reference 2. 


SYMBOLS 


a vector distance from point P in space to vortex element 

(fig. 1), ft 

A = x' cos t + y' sin i|r 

B = z' cos X - x' sin X 

C = (y' cos - x' sin ty) cos X - z' sin X sin i|/ 
ds vector length of vortex element, ft 

constant portion of wake vorticity, ft/sec 
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unit direction vectors along X-, Y-, and Z-axes, respectively 

running coordinate along edge of wake, ft 

arbitrary point in space 

vector induced velocity at P, ft/sec 

rotor radius, ft 

distance between P and edge of rotor disk at if (fig. 1), 

^R 2 + x 2 + y^ + z 2 - 2R(x cos \Jr + y sin if ) , ft 

dis tance betwee n P and center of rotor (fig. 1), 

Vx^ ~+~y2 + z^, ft 

vortex-element radius, measured parallel to tip-path plane 
from wake axis (fig. 1(b)), ft 

vector distance from origin to surface of cylindrical wake, ft 

induced velocity components directed parallel to X-, Y-, 
and Z-axes, respectively, ft/sec 

forward speed of rotor, ft /sec 

normal component of induced velocity at center of rotor with 
uniform radial and azimuthwise circulation, positive 
upward, ft /sec 

Cartesian coordinate system centered in rotor, X positive 
rearward, Y positive on advancing side of rotor, and 
Z positive upward 

Cartesian coordinates of point P (fig. 1), ft 
angle of attack of rotor tip-path plane, radians 

blade circulation, ft^/sec 

V sin a + Wa 
rotor inflow ratio, — 




rotor tip-speed ratio, V cos a 

'OR 



k 


X wake skew angle , angle between axis of tip-path plane and 

axis of wake, tan - -*- radians 

' A 

\|r azimuth angle , measured in direction of rotation from down- 

wind position (X-axis), radians 

H rotor rotational speed, radians /sec 

Primes on symbols denote nondimensionalisation with respect to R. 


THEORY 


Assumed Wake 

The shape of the wake assumed as a basis for the present calcula- 
tions is an extension of that used for previous rotary-wing induced- 
velocity calculations (refs. 2 and 3> for example). Its general char- 
acteristics are repeated herein for completeness. 

The circulation on the blades is considered uniform along the 
radius. (This assumption is in no way a restriction, since the results 
obtained under it can be converted to correspond with any arbitrary 
radial circulation distribution by the methods of ref. 1.) The free 
vortices leave the blade tip and are carried off with the speed and 
direction of the mean flow at the rotor. These vortices thus lie on 
the surface of an elliptic cylinder (fig. 1(a)). The vortices are 
assumed to be so closely spaced that the cylinder may be considered to 
be a sheet of continuous vorticity. This assumption restricts the 
analysis to obtaining only the time-averaged value of the induced veloc- 
ity (ref. 4). The axial component of vorticity is assumed to be negli- 
gible. This assumption is equivalent to stating that the rotor tip 
speed is infinite (ref. 5 ) y so that the circulation is proportional 
merely to the local disk loading. 

The essential addition in the present analysis is that the wake 
vorticity is allowed to vary in any arbitrary manner with the azimuth 
position. As a consequence of this azimuthwi r.e vorticity variation, 
the blades must also shed vorticity at their trailing edges. Figure 2 
illustrates this radial vorticity for a simple* case in which the circu- 
lation increases by an arbitrary constant over one sector of the rotor 
disk. In practice, however, the circulation changes are continuous, so 
that the wake cylinder will be completely filled with radial vorticity 
of a strength proportional to the derivative of the circulation. 
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The effect of the constant part of the blade circulation itself on 
the induced velocities is zero on a time-averaged basis. This result is 
not obtained, however, for the harmonic components of the blade circula- 
tion. Nevertheless, this contribution is ignored in the present analysis 
since it may be obtained relatively easily, and since it is zero in the 
only numerical example considered herein. 

The derivations which follow include all three components of the 
induced velocity. For convenience, the effects of the outer wake (the 
circumferential vorticity around the cylinder) and the inner wake (the 

L radial vorticity within the cylinder) are considered separately. 
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7 Induced Velocities of Outer Wake 

The induced velocities of the outer wake are found by integrating 
the Biot-Savart law over the entire wake. Thus, 


dq = A- ££ x a dL 
bn dL , _ , 3 
a 


(1) 


From figure l(a ) } "the following quantities may be determined by 
inspection: 

5 = i( R cos f + L sin X) + j(R sin \|r) + k(-L cos X) 
ds = i(-R sin \|r) + 3( R cos i|r) + k(0)J di|r 
a = i(R cos \|r + L sin X - x) + j(R sin t - y) + k(-L cos X - z) 

Substituting these values into equation (l) and integrating yields 


5-f 


at 


i J k 

-sin t|r cos ^ 0 

dr Kr cos j> + L sin X - x) (R si n f - y) (-L cos X - z) 
dL 


dL d+ 


[(R cos \|r + L sin X - x) 2 + (R sin + - y) 2 + (-L 


COS X - z 


>g 


3/2 


( 2 ) 



nent.- The normal (or k) component of Induced velocity 1b, from equation (2), 
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If ' now, g is expressed as where is the constant 

part of the vorticity, and where f(40 is a Fourier series in (nor- 

malized with respect to the constant term) , it is evident that the only 
term which can produce a normal induced velocity at the center of the 
rotor is the constant term. Thus, the induced velocity at the center of 
the rotor is 



so that finally 


r 2tT r 

2L = _L / f(^) U ~ ( x ' G0S ^ + y' sin t) + H c ' sin X cos ^jd^ 

W 0 ^ J o [ R c' + (cos i|r - x')sin X + z' cos X]R C * 


( 9 ) 


Longitudina l component .- The longitudinal (or i) component of the 
induced velocity is, from equation (2), 


r<=n r°° 

u = =2. f [ 

f dL 


(z + L cos X)cos dL d\|? 


Jo Jo [h 2 + x 2 + 


y 2 + z 2 - 2R(x cos V + y sin i|r) + 2L(z cos X - x sin X + R sin X cos i|r) + L 2 j 


3/2 


( 10 ) 


Tlie integration follows in precisely the same manner as that for 
the normal component. The final expression is found to be 


u 




(z 1 + R c r cos X)cos \|/ d^ 

[Rc* + (cos \|r - x' ) sin X + z' cos x] R c ' 


(ID 


Lateral component .- From equation (2), the lateral (or j) com. 
ponent is 


v r 


2rt r °° 

1 I fa r 


(z + L cos X)sln ^ dL d\|r 


Jo Jo [ r2 + x 2 + y 2 + z 2 - 2R(x cos i|r + y sin * ) + 2L(z cos X - x sin X + R sin X cos *) + L 2 ] 


3/2 


( 12 ) 


It is apparent that equation (12) differs from equation (10) only 
by a factor, tan X, which does not enter into the integration with 
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respect to L. Thus, the final expression for the lateral component 
of induced velocity may be written immediately, on comparison with 
equation (ll), as 

, r 2 * ( z ' + R ' cos X)sin ^ d+ 

= zi f(4r) - - ; 

w Q 2« J 0 R c • + (cos V - x')sin X + z' cos X|R C ' 


Induced Velocities of Inner Wake 


The induced velocities of the inner wake are found by integrating 
over the wake, where now, from figure l(b) , 

s = I(r cos V + L si n *) + j( r s ^- n + 4(-L cos 

ds = [i(cos V) + 3(sin *) + E(0)] dr 

a = i( r cos f + L sin X - x) + j(r sin i|r - y) + k(-L cos X - z) 

Substituting these expressions into equation (l) and integrating 
yields 




-at 


-R -<» 


0 -'0 


dtV 


I j * 

cob sin ^ ® 

(r cos >|r + L sin X - x) (r sin (-L co s x Z z ^- 


dL dr dV 


[(r cos ^ + L sin X - x ) 2 + (r sin \|f - y ) 2 + (-L cos X - z) 2 ] 


3/2 


(1M 


Normal component .- From equation (l4), the normal (or k) com- 
ponent of induced velocity is 


■Iff 


JLfEVx 

dt VdL/ 


sin ~ y cos >|f — L s: n X sin \Jf)dL dr d\|< 


Jo J 0 Jo [r 2 + x 2 + y 2 + z 2 - 2 r(x cos * + y sin *) + 2 L(z cs X - x sin X + r sin X cos *) + 


3/2 


(15) 


-q VO ->3 t" 1 



The Integration of equation (15) with respect to L is similar to the same integration for the outer wake {eq. and may be accomplished with the aid of items 162 and 170 of ref- 
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In certain cases, the argument of the inverse tangent term becomes 
indeterminate. The limiting value, however, may be found by the repeated 
application of L’Hopital's rule. Limiting values for a few of the more 
important cases are given in the following table: 


Location of P 

Azimuth angle, \|r, 
radians 

Arc tangent, 
radians 



, it 

n 

tsi 

n 

o 

0 

4 - _ 

2 




it 

0 


it 

X - 2 - 

x = z = 0 

2 

2 




3* 

2 

0 


Longitudinal component .- The longitudinal (or i) component of 
induced velocity is, from equation (l4). 


-i 

Ait 


jLflUVz + L cos X)sin if dL dr d| 
dt\dL J 


! 2 2 2 
l 0 + * + y 


2r(x cos if + y sin ijr) + 2L(z cos X - x sin X + r sin X cos \|/) + L 2 ] 




( 23 ) 


Equation (23) may be integrated with respect to L to yield 


c 2 * r 
-i I 


|^(^)[ Bin ♦][* + cos X^r 2 + x 2 + y2 + z 2 - 2r(x cos + y sin \|i)]dr 
' Jo Jo [)/r 2 + x 2 + y 2 + z 2 - 2r(x cos t|r + y sin \|i) + z cos X - x sin X + r sin X cos +] 'jr 2 +■ x 2 + y 2 + z 2 - 2r(x cos ♦ + y sin ^r) 

(2U) 


Integrating with respect to r, by means of equations (A5) 
and (Al4) of the appendix, results in 






Lateral component . - From equation (l4), the lateral (or J) component is 
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L 

7 

9 

7 
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is ma de in this paper to provide general flow charts for the various 
harmonics of vorticity such as has been done previously for the case of 
constant vorticity (for example, ref. 3). A few calculations for the 
normal component of induced velocity on the lateral axis of the roxor 
have been made for a unit sin \| t vorticity. These cases are discussed 
in the following sections. 


Distribution of Induced Velocity on Lateral Axis 

The calculated normal induced velocities are presented in figure 3 
for the case of the lateral axis and a unit sin V vorticity. The skew 

angle is tan”^" 4. The numerical results indicate that the induced veloc- 
ities are uniform on either side of the axis and discontinuous at the cen- 
ter of the rotor. This distribution corresponds to that of the assumed 
sin vorticity distribution, which is also uniform on either side of 
the axis but discontinuous at the center of the rotor. This close cor- 
relation between vorticity (or disk load) distribution and the induced 
velocities on the lateral axis was also noted in the experimental meas- 
urements of reference 1. 

The contribution of the inner and outer wakes individually is also 
shown in figure 3. The contribution of each part of the wake is a linear 
function of radial position - the outer wake portion increasing, and the 
inner wake portion decreasing with radius. Except for positions close 
to the edge of the rotor disk (y ' > O. 85 ), the inner wake produces the 
largest contribution to the total induced velocity for sin if vorticity. 

Reference 2 treats the induced velocities due to sin vorticity 
by means of a cruder approximation to the rotor wake. It is interesting 
to note that the actual induced velocity distribution (ref. 5) on the 
lateral axis, as obtained under the assumptions of reference 2, should 
be essentially the same, except in magnitude, as that obtained in the 
present analysis. The simple analysis of reference 2 cannot, however, 
be extended to locations other than the lateral axis. 


Effect of Skew Angle 

Since the absolute value of the induced velocity for sin ^ vortic 
ity is constant on the lateral axis, the variation with skew angle can be 
illustrated very compactly as in figure 4. Here the uniform value of 
induced velocity is shown as a function of tan X for a range of s ew 
angles encompassing the major portion of the aelicopter f g range, 
is evident that the effect of sin if vorticity increases somewhat with 
skew angle. The corresponding induced velocity from reference 2 for unit 

sin V vorticity is — = 2.0, irrespective of skew angle. Thus, even 

w 0 


— vo — J 
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though a strict interpretation of the procedures of reference 2 leads 

con^LST ^ape of the velocity distribution, it also leads to a 
siderable overestimate of the magnitude of the induced velocities. 


CONCLUDING REMARKS 


m «.h-i EqUat i° nG ^ WhlCh ° an bS lnte S rated by use of high-speed computing 
machines, have been developed for all three components of induced veloc- 

y at an arbitrary point near the rotor and for an arbitrary harmonic 
, h< 3 vor h lc ity distribution. Results of this investigation of the 
uced velocities near a lifting rotor with nonuniform azimuthwise 
/or lcity distribution show that for a vorticity distribution which 
varies as the sine of the azimuth angle, the normal component of induced 

rio shn y + a +h nS Slbher f de of the lateral ^is is uniform and antisymmet- 
nc about the center of the rotor. J 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., March 10, i960. 
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APPENDIX 

EVALUATION OF CERTAIN INTEGRALS OCCURRING IN THE ANALYSIS 


The analysis of the inner wake involves integrals of the form, 


dx 


X + (cx + d)\/X 


and 


I 


dx 


VX + cx + d 


where 

X = a + bx + x^ 

and where x is an arbitrary variable and a, b, c, and d are 
arbitrary constants . 

Both forms may be evaluated by use of the substitution (item 235 > 
ref. 6), 

u = N/X - x (Al) 

The first form is the simpler of the two and will be evaluated 
first. The indicated substitution (eq. (Al)) reduces this form to 


L 

7 

9 

7 


/ 


x + 


dx 

(cx + d )\jx 



2 du 

- l)u^ + (b - 2d)u t- (db - a - ac) 


(A2) 


The right-hand side of equation (A2) may oe integrated by item 67 
of reference 6 to yield 


P <ix _ 4 t -1 

2(c - l)u + b - 2d 

J X + (cx + d )]/x *pL 

ft J 


(A3) 


(A4) 


where 


q = 4db(c - 1) + 4a(l - c^) - (b - 2d)^ 
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The final expression is obtained by returning equation (A3) to the 
original variable x which results in 


dx 


X + (cx + d)\[x 


— tan"^- 

n 


2(c - 1)(\/X - x) + b - 

ft 


2d 


(A5) 


where q is as given previously in equation (A4). 

Evaluation of the second form is a more lengthy task but proceeds 
in the same manner. Making the same substitution (eq. (AX) ) as before 
yields 


dx 


2(bu - a - u^)du 


\/X + cx + d ( b - 2u)[u2( c - 1) + u(b - 2d) + (db - a - ac)] 


(a6) 


Now separate the right-hand side of equation (A 6) into partial 
fractions to obtain 


/ 


dx 

{x + cx + d 


i C 2 du + 2c r u du 

1 + c J b - 2u 1 + c J [u 2 (c - 1) + u(b - 2d) + (db - a - ac)J 


_ 2d r du 

1 + c J [u 2 (c - 1) + u(b - 2d) + (db - a - ac)] ^ 

' The first term may be integrated by inspection: the second term, 
by item 72 of reference 6 to give 

I W + » + d = IT7 l0ge(b ' 2u) + ^2TT log e[ u2 ( c - x ) + u ( b - 2d) + (db - a - ac)] 


be 


2d 



du 

- 1) + u(b - 2d) + (db - a 


ac)] 


c 2 . 1 


(A8) 
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The remaining integral may be evaluated by means of item 67 of 
reference 6 to yield 


P dx 

J Jx + cx + 


-1 


d 1 + c 


log P (b - 2u) 


c 4 * - 1 


log e ju^-(c - 1) + u(b - 2d) + (db - a - ac) 


5 2 - 1 


1- tan' 1 

2(c - l)u + b - 2d 

n 

£ J 


(A9) 


where q is defined in equation (A4) . 

Returning equation (A9) to the original variable x results in 


P = — = log e (b + 2x - 2\/x) 

j^C+cx+d c 2 - 1 


;2 . 1 


logp 


(7x - x) 2 (c - 1) + (\/x - x) (b - 2d) + (db - a - 


bc - 2d 2. tan' 1 


j~ 2(c - 1) ( l/x - t) + b - 


=2-1 ^ 


2d 




(A10) 


Note that 

(\/x - x) 2 (c - 1) + - x)(b - 2d) + (db - a - ac) = (b + 2x - 2 \Zx)(\fc + cx + d) 

so that the logarithmic terms of equation (A10) may he combined to 
yield 


P dx 

J + cx + 


d - 1 


log e (\/x 


+ cx + d) + — log e (b + 2x - 

c*~ - 1 


be - 2d 2 , -1 

tan 


- 1 ]/q 


2(c - 1 ){\fx - x) + b - 2d 


(All) 
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The presence of the logarithmic terms in equation (All) occasions 


imaginary values for the integral if either (tfx + cx + dj or 
(b + 2x - 2 'jx ) is negative. Provided that these terns are monotonically 


+ cx + d 


) 


negative, this result may be avoided as follows: 
Rewrite equation (A7) as 


r dx 

1 

f 2 du 

2c [ 

J VT + cx + d 

1 + c 

j 2u - b 

1 + c j 


h 


u du 


[u 2 (l - c) + u(2d - b) - (db 


2d 

1 + c 


du 


[u 2 (c - 1) + u(b - 2d) + (db - a - ac)] (A12) 

Integration, as in equation (a 8), of the first two terms results in 

f - j= — = ^ — lo K e ^ 2u - b ) + -p" 2 log [u 2 (l - c) + u(2d - b) - (db - a - ac)l 

J U + cx + d 1 + c - 1 L J 


be - 2d 
c 2 - 1 



du 

1) + u(b - 2d) + (db - a - ac)j 


(A13) 


Comparison of equations (A8) and (A13) indicates that the only 
difference is a reversal in sign of the arguments of the logarithmic 
terms . Therefore, as long as these arguments are monotonically either 
positive or negative, the final expression may be rewritten as 


dx 


Vx + cx + d 


c 2 _ lo 6e |Vx + cx + d 


+ — lo s e | b + 2x 

c 2 - 1 e| 



where , again. 


5!S_r^i 2. ^-1 

c 2 - 1 ^q 


2(c - 1) (& - x) + b - 


2d 


& 


(A14) 


q = 4db(c - 1) + 4a(l - c 2 ) - (b - 2d) 2 
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Figure 4.- Absolute value of normal component of induced velocity on either side of lateral (Y) 

axis for unit sin vorticity. 




